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Summary 

The  great  success  of  the  two  S1R-C/X-SAR  missions 
and  the  lack  of  global  topographic  data  are  the  reasons 
at  the  basis  of  third  Shuttle  Radar  Topographic  Mission 
(SRTM),  The  mission  is  a joint  program  among  NASA 
NIMA,  DLR  and  the  Italian  Space  Agency.  The  SRTM 
mission  aims  is  to  create  a tri-dimensional  topographic 
map  of  the  Earth's  surface  30  times  as  precise  as  the 
best  global  maps  in  use  today. 

The  sophisticated  mission  technology  is  able  to  produce 
topographic  maps  with  a level  of  detail  never  achieved 
before  and  a very  good  height  precision. 

The  technique  used,  based  on  microwave  instruments 
with  a synthetic  aperture  radar,  is  a totally  new  concept 
in  geodesy,  the  science  that  studies  the  morphology  of 
planet  Earth  and  measures  its  characteristics. 

The  mission  has  been  also  a great  opportunity  to  test  in 
an  operational  contest  a innovative  one  bit  processing 
chain  operating  in  the  pure  time  domain.  The  processor 
focuses  the  X-SAR  signum  coded  raw  data  applying  a 
quasi-phase  preserving  algorithm  The  equivalent 
computing  power  is  about  8,000  billions  operations  per 
second. 

The  processor  architecture  is  such  to  comply  with  the 
continuous  data  flow  coming  from  the  instrument 
delivering  a striped  image. 

Alenia  Spazio.  beside  the  responsibility  with  Dornier  of 
the  design  and  realisation  of  the  X-band  radar,  also  in 
this  case  under  ASI  contract,  has  designed  and  realised 
a complete  ground  station  organised  as  a X-band  data 
Processing  Chain  installed  and  operated  at  JPL  site. 
This  includes  a Real-Time  SAR  processor,  a 
comprehensive  data  management  system  and  finally  a 
topographic  post-processing  station. 

Astonishing  results  have  been  produced  both  as  single- 
bit focused  SAR  images  and  digital  elevation  models, 
also  thanks  to  the  advanced  capabilities  developed  in 
last  years  at  University  Federico  II  of  Naples  one  of  the 
more  active  research  institutions  in  this  field. 


Background 

X-SAR,  X band  Synthetic  Aperture  Radar,  was 
developed  by  Alenia  Aerospazio  for  Italian  Space 
Agency  (ASI),  as  part  of  the  joint  Italian,  Gentian  and 
US  space  Shuttle  Radar  Laboratory  (SRL)  program. 

It  was  designed  for  integrated  operation  with  the  SIR-C 
Synthetic  Apenure  Radar  developed  by  JPL  as  part  of 
the  same  mission. 

The  instrument  was  successfully  qualified  and  validated 
through  two  space  shuttle  missions  in  1994,  were  its 
challenging  imaging  performance  at  X band  were 
proven. 

It  was  the  first  X-band  SAR  developed  for  spacebome 
applications,  and  provided  a significant  contribution  to 
improve  the  knowledge  of  multi-frequency  SAR 
imaging  in  co-operation  with  the  L/C  band  of  SIR-C 
instrument. 


New  Features 

X-IFSAR  is  a radar  interferometer  based  on  a couple  of 
Synthetic  Aperture  instruments,  that  share  the 
transmitting  antenna  while  hold  separate  the  receiving 
channels  including  antennas,  RF  sections  and  digital 
data  handling  subsystems. 

The  X-Band  acquisition  is  guaranteed  providing  all 
weather  operation  and  night/day  imaging  and 
topographic  mapping  of  about  35%  of  the  Earth  surface. 

The  Italian  Processing  Chain  operated  in  connection 
with  X-IFSAR  in  Pasadena  (CA)  at  the  Jet  Propulsion 
Laboratory  NASA  centre. 

The  direct  down-linked  data  w'as  available  for 
processing  and  delivering  first  absolute  DEM  products 
together  with  auxiliary  functions  aimed  to  instrument 
in-flight  calibration  and  performance  monitoring. 

The  X-IFSAR  Processing  Chain  employees  a 
revolutionary  real-time  parallel  processor  based  on  a 
proprietary'  ASIC  solution.  Each  of  the  128  ASICs, 
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developed  using  VHDL  coding  and  Field  Programmable 
Gate  Array  (FPGA)  technology,  endorses  the  Single  Bit 
SAR  coding  and  focusing  algorithm,  first  conceived  and 
analysed  by  the  University  of  Naples  "Federico  II”, 
Department  of  Electrical  Engineering. 

The  Signmn  Code  (SC)  algorithm,  together  with  state- 
of-the-art  FPGA  chip  architecture,  made  the  real-time 
operation  possible,  achieving  full  geometric  resolution 
on  SAR  images  and  full  height  accuracy  for  DEM 
generation. 


Figure  l The  TD  Complex  Correlator  FPGA 


Applications 

X-IFSAR  is  derived  from  X-SAR  whose  primary 
application  has  been  Earth  surface  imaging  together 
with  the  study  of  the  moisture  connected  to  hydrology, 
geology,  ecology'  and  other  disciplines. 

The  specific  purpose  of  X-IFSAR  is  the  measurement  of 
ground  height  in  order  to  deliver  high  resolution,  high 
accuracy  topographic  maps  all  over  the  world.  The 
planning  of  the  11-davs  SRTM  mission  constrained  X- 
IFSAR  to  acquire  about  35  % of  the  Earth  surface 
within  +/-  57°  latitude  boundaries. 

The  Italian  Operational  Processing  Chain  has  been  used 
within  the  mission-time  frame  for  in-flight  calibration 
purposes  and  absolute  first  product  generation. 

After  SRTM  conclusion,  it  has  been  installed  al  "Centro 
di  Geodesia  Spazialc”  Matcra  - Italy  - where  all  the 
SRTM  data  set  has  been  available  starting  from 
November  1999.  There  the  Processing  Chain  has  been 
used  for  systematic  fast  product  generation. 

The  comparison  with  the  large  data  sets  acquired  during 
SRL  missions  in  1994  will  also  allow  the  observation  of 
midterm  geological  process,  soil  erosion  and  deserts 
progression  within  a five  years  period. 


The  huge  coverage  granted  during  SRTM  by  X-IFSAR 
will  allow  once  again  the  analysis  of  precipitation  over 
tropical  forest  and  ice  properties  and  their  variation  and, 
finally,  ocean  large-scale  dynamics 


Basic  SC  Arithmetic 

The  arithmetic  operation  that  plays  the  central  role  in 
the  novel  SC  data  representation  is  the  multiplication. 

Given  h and  r two  variables  defined  within  the  real  field 
number  (usually  represented  adopting  the  1EEE754 
standard  on  conventional  computer  FPUs),  two  new 
variables  can  be  associated  as  defined  here  below 


h->  h = 


h > 0 
h < 0 


and 


J+l  r>  0 
{-1  r<0 


The  new  variables,  generated  by  means  of  a hard 
quantization,  also  belong  to  the  real  number  field,  but 
they  can  assume  only  two  values. 

Now  the  new  variables  can  be  re-coded  using  a more 
appropriate  representation  that  accounts  for  true 
possible  states.  The  binary  representation  is  immediate. 


Let’s  //  and  R be  the  codes  corresponding  to  " and  r : 


h-*H 


fl  A = +1 
|0  A = -1 


and 


f -»•  R = 


r = +l 
r = — 1 


The  binary  coding  of  hard  quantized  variables  H and  R 
is  named  Signum  Coding  shortened  in  this  frame  with 
the  symbol  SC. 

Now  let's  examine  how  the  signum  coding  works  on 
arithmetic  multiplication 


47-3 


The  h by  r product  generates  the  results  s that  can  be 
represented  using  binary  coding  too. 

It  can  be  shown  that  the  value  of  .S',  the  SC  of.?,  obtained 
by  means  of  the  negation  of  an  exclusive  OR  between  // 
and  R.  is  exactly  equal  to  the  value  got  applying  SC  on 
the  hard  quantization  of  s-  h r . 

In  term  of  mathematics  this  property  can  be  written  as: 
s = h-r=*S  = R<2>H 

which  can  be  verified  case  by  case  looking  into  Table  1, 
where  S,  computed  in  the  two  equivalent  ways,  are 
systematically  compared. 


Table  1 - Real/Boolean  relations  for  multiplication 

If  the  True  State  Counter  (TSC),  or  “1”  occurrences 
counter,  is  defined  using  the  symbol 


The  time  domain  operation  used  for  SAR  focusing  is  the 
bi-dimensional  correlation,  performed  using  mono- 
dimensional correlation  followed  by  results  addition 
This  correlation  is  defined  in  the  complex  number  field 
and  can  be  decomposed  into  four  real  number  field 
correlations. 

Therefore  the  basic  function  is  the  real  correlation, 
which  is  susceptible  of  signum  coding  representation. 

Real  Correlation 

Let  r(*,*)  be  the  input  complex  data  matrix  built  of  the 
aligned  raw  data  radar  echoes,  and  h(**)  be  the 
complex  bi-dimensional  space-varying  holographic 
SAR  reference  function. 

The  derived  measurement  .?(*,*)  of  local  ground  whet 
reflectivity  can  be  computed  in  the  native  time/space 
domain  as  the  correlation  product : 

Nr  Na 

s(i,  y>  = z S r(‘  + pJ+< 7)  • hi.j  (/>.?)- 

p~\q-^\ 

The  space  varying  dependence  of  is  depicted 

by  the  i and  j indexes  The  hard  quantization  transforms 
the  expression  of  .?(/,  /)  into 

Nr  Na  A 

KiJ)  = XX Hi  + PJ  + ‘ 7) • h,j (p>(])  ■ 


\z 

H 

R®  II 

s | 

1 

1 

1 

+ 1 

1 

0 

0 

-1 

0 

1 

0 

-1 

0 

0 

1 

+1 

r ■ h 


+1 

+1 

+1 

+1 

-1 

-1 

-1 

+1 

-1 

-1 

-1 

+ 1 

l|{.}  counts  the  number  of  " Sj  = 1 

J Wj  cnt  occurrences  while  j varies 


the  sum  of  products  x,  can  be  written  using  the 
equivalent  SC  coding 


It  is  not  relevant  now  how  the  s(ij)  and  s(i,j) 
expressions  differ,  because  the  theoretical  SAR  focusing 
performance  is  not  a matter  of  interest  here 

If  the  equivalence  properly  of  the  sum  products  is 
applied,  the  s(',j)  can  be  immediately  evaluated 
employing  the  SC  representation  and  getting  the 
relationship 


= 2 ■ S’ 


J II  j.CNT 


- N 


where  N,  the  number  of  added  elements  x,,  acts  as  a 
bias.  N is  sometimes  called  “magic  number”. 


Signum  Code  SAR  Focusing  Algorithm 

Here  the  SAR  processing  algorithm  based  on  the  SC 
raw  data  representation  is  detailed  as  implemented 
within  the  X-SAR  Processing  Chain. 


.?(/,  j)  = 2Z  + P-  j + 9)  ® Hu  (Pi  <7)1^ 

-NpN 


<7 


The  bi-dimensional  count  operation  is  exploited  into  the 
sum  in  one  dimension  of  mono-dimensional  counts  in 
the  other  one. 


47-4 


Complex  Correlation 

The  complex  correlation  is  similar  to  the  real  one 
analysed  in  the  previous  section.  The  mathematical 
definition  of  it  is  exploited  by  the  following  equation 


Nr  Na 

s(ij)  = £ £ r(i  + p,  ;+<?)■  h' , ( p,  q) 


/'=«<7=1 


Now,  by  applying  the  discussed  niles  and  dividing  by  2 
(which  does  not  break  the  generality'  of  this  property), 
the  SC  representation  S of  s can  be  computed  as 


C-^-t 

2 


= pr 


®/c 


+\\aT  ®b" 


c"  = -a 


: 


where  h,  f(p.q) , r(/,y)  and  .s(/,  /)  are  complex 
numbers  instead  of  real  numbers. 

Note  that  the  reference  function  h,  , (p,  q)  shall  be 

conjugated  in  order  to  express  a full  significant 
correlation. 

The  complex  multiplication  between  two  complex 
numbers  can  be  reduced  into  a set  of  four  real 
multiplications  and  their  proper  combination  for 
generation  the  real  and  imaginary  part. 

Given  A and  B the  signum  code  representations  of  two 
complex  numbers  where,  ARe.  BRe,  Ata  e Bu"  stay  for  the 
four  components,  the  complex  multiplication  of  A by  B* 
(*  signifies  complex  conjugation  function). 


Note  that  the  imaginary  part  CRc  only  requires  the  bias 
compensation  due  to  the  presence  of  positive  counters. 
On  the  other  hand  the  biases  of  the  two  C11"  components 
compensate  each  other. 

For  the  purposes  of  the  SAR  focusing  application,  the  by 
2 division  does  not  break  the  validity  of  the  property. 

Non  space-varying  SAR  focusing 

The  ideal  SAR  focusing  scheme  has  been  employed  that 
fully  respects  the  intrinsic  space-variance  of  SAR 
acquisition  system.  In  fact  the  the  reference  correlation 
function  )i  ,{p,q)  depends  on  the  raw  data  sample 

position  through  i and  j. 


C = A ■ B~ 


is  equivalent  to  the  combination  of  six  real  field 
operations 


C-> 


| rRe  -RRe 

]clni  =-aKc  -b"" 


f ,4lm  • Blm 
+ A,m-BRt 


The  implementation  of  full  space-variant  filter  is  not 
realistic  and  would  lead  to  an  extremely  complex  design 
resulting  not  compatible  with  real-time  architecture  and 
run. 

The  relationship  that  determines  the  translation  from 
raw  data  to  focused  image  shall  be  therefore  simplified 
in  some  way. 

Keeping  the  same  meanings  for  >■(*,*)  and  .v(*,*)  an 
alternate  relation  is  derived  that  expresses  s(i,  j)  as 


This  relates  to  the  complex  product  only.  If  the  complex 
addition  is  introduced,  the  sum  of  complex  products 


Nr  Na 


j ) = I X r(i  + p.  j + q)  ■ hK  ( p,  q) 

p=\  q=\ 


can  be  translated  in  the  real  field  according  to  the 
scheme  here  below 


=-2>f  - i>t+ 


where  hK(p,q)  is  frozen  in  connection  with  a limited 
number  of  K values  (K  from  1 to  5 are  reasonable)  in 
such  a way  to  assume  a piecewise  constant  reference 
function  The  validity'  region  of  each  filter  make  the 
processing  act  as  ’unfocused’. 

The  requirement  baseline  requests  for  one  single 
correlation  filter  all  over  the  swath  that  shall  be 
matched  in  correspondence  of  the  swath  centre.  The 
azimuth  focusing  degradation  implied  by  that  choices 
directly  proportional  to  the  swath  width  where  the 
single  filter  operates  a space-invariant  correlation 
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Figure  2 Time  Domain  Complex  Correlator  Architecture 


Fast  Filter  Evaluation 

If  rcs  is  (he  centre  swath  range,  that  is  the  distance 
between  the  antenna  phase  center  and  the  swatli  center, 
the  Doppler  rate  pA  induced  by  the  effective  relative 
motion  is  given  by: 


where  the  negative  sign  is  due  to  the  physical  change  of 
radar-to-target  distance  during  SAR  motion. 

The  azimuth  quadratic  phase  excess  d> ^ (f),  as  function 

of  flight  time  t computed  at  swath  center,  is  then 
defined  as 

= o)2  Wlth  |^'o|<Y 


where  f0  is  the  flight  time  over  a specified  target  placed 
at  swath  center  and  TA  is  the  coherent  integration  time. 
i is  the  flight  time. 

This  is  an  exact  expression  at  swath  centre  only  and 
should  be  corrected  for  any  other  range  position  The 
assumption  of  using  a single  reference  filter  consists  in 
extending  the  validity  of  this  quadratic  phase  shape  to 
all  ranges  within  the  swath. 

The  total  azimuth  phase  excess  shall  include  also  the 
linear  term  induced  by  the  Doppler  centroid,  physically 
due  to  a slight  skew  in  acquisition  geometry.  If /Dc  is  the 
Doppler  centroid  frequency  the  azimuth  linear  phase 
term  is  given  by 

(1>a/(0=2^-/Cc(f-f0)  wilh  Mo|<Y 


Tlie  range  quadratic  phase  is  instead  determined  by  the 
chirp  pulse  bandwidth  and  duration  If  Tp  is  the  pulse 
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duration  and  Bp  its  bandwidth,  the  chirp  rate  is 
simply  given  by  the  ratio 


Mn  =y~ 

1 P 

then  the  range  quadratic  phase  d>/((r).  as  function  of 
echo  fast  time  x.  can  be  computed  as 

T 

®„(r)=  w"h  |r-ro|<y 


where  xo  is  the  centre  of  the  chirp  pulse. 

The  correlation  filter  hK(p,q)  (this  is  not  the  signum 
coded  filter  but  the  full  dynamic  one)  is  defined  over  the 
2D  domain  determined  by  the  range  axis  (p  index 
samples)  and  the  azimuth  axis  (q  index  samples). 

The  relation  between  range  index  p and  echo  fast  time  x 
is  defined  by  the  sampling  frequency  fs 


while  the  relation  between  the  azimuth  index  q and  the 
flight  time  t is  determined  by  the  pulse  repetition 
frequency  (PRF)  as 


l — — “ + l fjpi' 

1 PRF 


where  twa?  is  arbitrary  and  can  be  ignored  without 
lacking  of  generality. 

Assuming  the  following  ranges  for  p and  q: 

p = 1,2, , Nr 

q = 1,  2 Na 

the  time  centres  to  and  t0  are  assumed  in 
correspondence  of 

1+  N 

t„  = ~ 

2 • PRF 


and 


l + N 


so  that  the  total  phase  ^(t,r)  over  the  filler,  that  can  be 
computed  in  the  2D  continuous  domain  as 

<Kr-0  = ^.,(0  + ^/(')  + «Ur) 

T T 

with  |r  - r J < — and  /-/„<  — 

1 1 2 ' 01  2 

is  re-written  according  the  expression 

7>{p*cl)  = h<,{(l)  + 7 Ai  (q)  + 7k  (p) 

with  p = l,..,Nr  and  q~l,..,Na 
where 


From  the  total  phase  <j)(p,q)  the  real  and  imaginary 
parts  of  the  correlation  filter  hK(p,q)  can  be  computed 

as: 

h*'-(p,q)  = cos[0(p,  <7)] 

h'l?{.p,q)=seif${p,q'h 

which  are  the  full  dynamic  versions  of  the  correlation 
filter  required  by  Signum  Code  SAR  focusing 
algorithm. 

The  Signum  Code  version  can  be  immediately  derived. 
Anyway,  without  computing  the  trascendent  functions. 
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it  is  possible  to  generate  the  signtim  of  cosine  and  the 
signum  of  sine  using  the  simple  verification  rule  based 
on  the  phase  normalization 

A normalized  phase  4'  is  associated  to  the  true  phase  0 
by  means  of: 


4'  = mod 


l<t> 


.4 


that  furnishes  4/  in  the  range  |0,4[.  Then  the  sine 
signum  Ssulc  and  cosine  signum  Scosm,  are  immediately 
obtained  looking  for  the  integer  value  of  4J  achieving  at 
the  end: 


Z = int(Y)  -> 


^cosine  ^ ^ sine  ^ ^ 

cosine  ^ sine  ^ ^ 

^ cosine  ~ ^ ^ sine  ~ ^ 

S ,,,=+!  .S'  =-l 


if  Z = 0 
if  Z = 1 
if  Z = 2 
if  2 = 3 


System  Operation 

The  system  operation  was  based  on  joint  operation  of  an 
On-Board  Segment,  including  the  SAR  sensor,  with  the 
Ground  Segment,  including  the  necessary  equipment  for 
mission  planning,  control  and  data  image  and  DEM 
processing. 

The  On-Board  segment  includes  the  two  X-band 
antennas  (12  in  and  6 m azimuth  apertures),  integrated 
with  the  C-band  antennas,  the  RF  and  digital 
subsystems  for  interferometric  radar  operations.  Digital 
data  was  formatted  and  stored  on-board  and  a fraction 
(~2%)  was  down-linked  to  the  Jet  Propulsion 
Laboratory. 

The  Italian  Processing  Chain  for  X-1FSAR  was 
supported  by  the  mission  and  sensors  planning 
operators.  Most  of  instrument  parameters,  significant 
for  final  image  and  DEM  performance,  was  monitored 
in  conjunction  with  the  German  and  US  Ground 
Segments. 

The  X-IFSAR  Processing  Chain  supplied  the 
acquisition  and  monitoring  of  housekeeping  data 
received  from  the  sensor  via  the  high-rate  bit- 
multiplexed  channels,  as  well  as  the  real-time  focusing 
processing  of  received  scientific  data  through  SAR 
images  generation.  All  bit-stream  synchronisation  was 
performed  using  latest  hardware  and  software  available 
technologies. 


Figure  3 Ch’erall  SAR  focusing  SC  algorithm  block  diagram  implementation 


47-8 


Key  Technologies 

The  emphasis  of  X-SAR  development  was  put  on 
application  of  the  latest  available  space  qualified 
technologies. 

Now  the  technologies  employed  within  the  advanced  X- 
IFSAR  Italian  Processing  Chain  are  the  key  through  the 
incoming  Italian  Space  Programs,  such  as  COSMO  and 
SAR2000,  that  will  benefit  of  the  drawback  of  the 
developments  and  results  achieved  in  the  frame  of 
SRTM  program. 

The  present  SAR  processor  is  a machine  that  sizes  as 
large  as  one  VME  crate  cabinet  (8  boards)  that,  thanks 
to  the  acquired  know-how,  will  be  miniaturised  down  to 


a single  double-Europe  card  and  few  tens  of  Watt  power 
supply. 

These  feature  and  specifications  in  terms  of  envelope, 
mass  and  power  enables  the  on-board  implementation  of 
the  Siginim  Code  SAR  processing  technique  which  is 
seen  now  as  the  switch-on  gate  toward  new  applications 
of  spaceborne  SAR  instruments  addressed  to  a larger 
base  of  users. 

This  is  the  key  relevance  of  the  Italian  participation  to 
SRTM  where  SRL  experience  and  academic  invention 
has  been  exploited. 


Figure  4 Owrall  Processing  Chain  Architecture  as  implemented  at  JPL  site 


Data  Management  Service 
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Here  are  listed  main  specifications  of  the  X-SAR 
Processing  Chain  components. 


SAR  Processor 


implementation 

Full-custom 

SAR  Processor  Algorithm 

Time-Domain  Signum  Code 

Operation  Capability 

Real-Time 

Architecture 

SIMD 

CPU  Technology 

FPGA / BGA 

Number  of  CPUs 

128 

Clock  Frequency 

30  MHz 

Equivalent  Computing  Power 

- 6,500,000  Mops 

Master  Controller 

Motorola  68060 

RT  Operating  System 

OS9k 

Input  l/F:  raw  data 

IEEE  1284 

auxiliary 

RS232C 

command 

VME-PCI  bridge 

Output  Interface 

Fast  Ethernet 

Local  Storage 

63  GByte 

Recording  Autonomy 

35  minutes 

Host  Computer 

SUN  E450  WS 

Disk  Storage 

117  GByte 

RAM  Memory 

1 GByte 

CPUs 

4 x 400  MHz  - UltraSPARC 

I/O  Interfaces 

Fast  Ethernet 

Tasks  on  X-SAR  raw  data 

Synchro-Deformat  - DQA/PE 

Service  Capability 

60  minute/day 

Network  Services 

FTP  Server  & WEB  Server 

Topographic  Processor 


Host  Platoform 

SUN  E450  WS 

Disk  Storage 

33  GByte 

RAM  Memory 

512  MByte 

CPUs 

4 x 300  MHz  - UltraSPARC 

I/O  Interface 

Fast  Ethernet 

Products 

Fringes  & DEM 

Processing  Capability 

1 DEM  / 2 hours 

Algorithm  & S/W  Source 

University  of  Naples 

Chain  Controller 


Host  Computer 

IBM-PC 

Platform  Architecture 

Intel  Pll  @450  MHz 

Operating  System 

NT  4.0 

Control  Environment 

LabView  5TM 

Local  Control 

full 

Remote  Control 

operations 

Thermal  Control 

on  128  CPUs 

Fan  Control 

VME  fans 

I/O  Control  Interface 

VME 

RT  Image  Data  Display 

Viper550-TNT 

On-Screen  Resolutions 

800x600 

1280x1024 

1600x1200 

VCR 

PAL  (VGA) 

System  Pnnter 

InkJetTM 

Photographic  Support 

NP1600 
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